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Abstract: In this paper we propose a semi-Markov modulated model of interest 
rates. We assume that the switching process is a semi-Markov process with finite state 
space E and the modulated process is a diffusive process. We derive recursive equations 
Oh \ for the higher order moments of the discount factor and we describe a Monte Carlo al- 

gorithm to execute simulations. The results are specialized to classical models as those 
U—\ ■ by Vasicek, Hull and White and CIR with a semi-Markov modulation. 
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^ ■ 1 Introduction. 

■ The literature on interest rate models is ample and mainly concerns models based on 

CN ■ short rate dynamics and models of forward rate, see e.g. |Bj6rk(2004)]. The advantages 

and the drawbacks of short rate models are well known in literature and even nowadays 
they receive continue attention by researcher and practitioners. 
The fundamental paper by [Vasicek( 1997)1 is the 

pioneering contribution of short 
rate models. Different papers presented extensions and alternatives to the Vasicek 
model. One interesting approach in term structure modelling is that of |Duffie and Kan(1996)| 
where it is assumed that the short rate process is a function of a state process. Inter- 
esting developments were subsequently obtained by |Mamon(2002)| who characterized 
the term structure of a Markov interest rate model when the interest rate process is 
assumed to be a function of a continuous time non-homogeneous Markov chain. By 
using the forward measure in [Mamon(2004)1 it was shown how to price term structure 
derivative products. 

The present paper adopts the assumption that the short rate process is a diffusive 
process modulated by a continuous time semi-Markov process. The process of interest 
rates is defined by using the theory of semi-Markov reward processes with initial back- 
ward times as developed by [S tenberg et al.(2007)| and here opportunely extended to 
consider the case of stochastic permanence rewards. In this way we provide a general 
model of evolution which is able to reproduce a great variety of evolutions of the interest 
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rate. The results include renewal type equations for the higher order moments of the 
zero coupon bond process and for the covariance function of the force of interest. No- 
tice that, in the paper by |Hunt and Devolder(20irj] , a discrete time regime switching 
binomial-like model of the term structure, where the regime switches are governed by 
a discrete time semi-Markov process, is presented. There, it is possible to find a clear 
description of the advantages of the semi-Markov approach. 

It should be highlighted that our paper doesn't discuss the problem of transformation 
between the real-world (physical) measure and the pricing measure and for this reason 
the proposed model could be adopted not for pricing interest derivative products but for 
all applications that require the real- world evolution of rates, see |Rebonato et al.(2005)| 
for a complete list. Finally notice that particular cases of our model were shown partic- 
ularly adapted to solve problems with applications to insurance, see [Norberg(1995)] . 

The paper is organized as follows. The next section present a short description of 
semi-Markov processes with the introduction of the notation adopted in the analysis. 
Section 3 presents the stochastic models of the short interest rates. In this section we 
derive the main results concerning the equations for the higher order moments of the 
zero coupon bond process and for the covariance function of the force of interest. Sec- 
tion 4 describes the results specialized to popular diffusive models of short rates and 
present a Monte Carlo algorithm able to generate the synthetic data of the model. 

Let us consider two sequences of random variables defined on a complete, filtered 
probability space (f2, J 7 , J-j, P): 

(i) J n : Q I = {1, 2, m}, n G N representing the state at the n-th transition; 

(ii) T n : Q — > R + representing the time of the n-th transition. 

We suppose that (J n ,T n ) is a Markov renewal process of kernel Q = [Qij{t)\: 

FVn+i =j,T n+ i -T n < t\a(J a ,T a ),J n = i,0 < a < n] 
= P{J n +i = j,T n+1 -T n < t\J n = i} =: Qtj(t). 

We also know that: 

Pij := P{J n +i = j\ J n = i} = hm Qij(t). 

t— »oo 

P = [pij] is the transition matrix of the embedded Markov chain {J n }nelN- 

We assume that t — » Q(i), t G M+ is absolutely continuous with respect to the Lebesgue 

measure on R + . 

Denote the probability the process will leave state i within time t by: 

Hi(t) := P{T n+1 -T n <t\J n = i} = ^2 Qik(t)- 

kei 

It is possible to define the distribution function of the waiting time in each state i, 
given that the state successively occupied is known: 

Gij(t) = F{T n+1 -T n < t\J n = i, J n+1 = j} (1.1) 

The related probabilities can be obtained by means of the following formula: 



Gijit) 



Qij{t)/Pij, if Pij ¥= o, 
1, if Pij = 0, 
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Denote by N(t) := sup{n : T n < t}, the semi-Markov process Z = (Z(t), t € R.+ ) 
is defined as Z(t) := Jjv(t)- 

The transition probabilities of Z(t) are defined in the following way: 

</Hj(t) :=F{Z(t)=j\Z(0)=i,T N{0) =0}. (1.2) 

They are obtained by solving the evolution equation: 

^•(t) = <%(1 - Hi{t)) + Vf Q ik (e)<f> kj (t - 9)d9 (1.3) 

fee/ ^ 

where <5jj represents the Kronecker symbol. 

The part 6ij(l—Hi(t)) gives the probability that the system does not have transitions 
up to time t given that it starts in state i at time 0. 

The term Ylkel Jo Qik^^kjif, — 6)d6 considers the permanence of the system in 
state i up to the time 9 where a transition in state k occurs. After the transition, the 
system will move to state j following one of all the possible trajectories going from state 
k to state j in the remaining time t — 0. All possible states k and times 9 are considered 
by the summation and the integration. 

Semi-Markov processes are very convenient to describe phenomena which display 
a duration effect. The duration effect affirms that the time system spent in a state 
influence its transition probabilities. One way to detect and quantify this effect with 
semi-Markov processes is by using backward and forward recurrence time processes. 
Recurrence processes were analyzed in |Janssen and Manca(2007)| and more recently in 



[D'Amico et al.(2009)| and [D'Amico et al.(2011) 



For our purposes we consider only the effect of the initial backward value. 
Given (J n ,T n ), we define the following recurrence process: 

B(t):=t-T N{t) . (1.4) 

B(t) is called the backward time (or age) process. 

The transition probabilities of a semi-Markov process change in function of the values 
of the backward time. In fact the conditional waiting times distribution functions (II. ip 
can be of any type and then, also no memoryless distributions can be used. In this 
case the time length spent in the starting state (initial backward value) changes the 
transition probabilities. 

The probability with starting backward b (frij(u; t) := P[Z(t) = j\Z(0) = i, B(0) = u] 
denotes the probability of being in state j after t periods given that at present the process 
is in state i and it entered into this state with the last transition u periods before. It 
results that 



<fa-M) = % 1- Hi {u) + ^i l- Hi (u) Mt - d)dd - (L5) 



If -B(O) = 0, then u = and equation f| 1 . 5 1) degenerates in (|1.3p . 
The probability with starting backward is a very special case of the transition prob- 
abilities studied in |D'Amico et al.(2009)| and |D'Amico et al.(2011)| . 
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2 The semi-Markov modulated model 

In this section we define a semi-Markov modulated model of interest rates and we 
assess its probabilistic behaviour. 

In the following we assume that the force of interest at any time t is a stochastic 
process of diffusive type whose parameters depend on the state of the semi-Markov 
process, on the backward recurrence time process and on the initial value of the force 
of interest ro- To be more precise we assume that within two transition times T n _i 
and T n of the semi-Markov process the evolution of the force of interest follows the 
dynamic of a diffusive process whose parameters depends on the state J„_i of the semi- 
Markov process. The dynamic of the force of interest, between two consecutive renewal 
moments, will be indicated as 

dr{t) = bi(r(t),t)dt + ai(r(t),t)dW(t), r(0) = r (2.1) 

where it is supposed that To = 0, T\ > t and Jq = i. The process W(t) is a brownian 
motion with respect its own filtration , and we denote by b and a the drift and the 
diffusion coefficients, respectively. The solution of the stochastic differential equation 
(|2,ip will be denoted by rj iro (s), for s € [Tq,T\). In T\ the semi-Markov process transits 
to another state, say j, the force of interest in the time interval between T\ and the 
next transition will evolve according to 

dr(t) = bjir^^dt + a j (r(t),t)dW(t), r(T x ) = r iinj (Ti) (2.2) 

for t E [Tij T2), and so on, then the resulting force of interest will be a continuos process. 
We would like to stress that for any s £ [Tq,Ti), and in general between any couple of 
transition times, the solution rj jro (s) is obtained with the standard methods of stochastic 
calculus. 

Due to the fact that the force of interest depends on the modulating process, we need 
to describe the force of interest process at any time s given the information available at 
the present time (time zero as long as we work with an homogenous time model) given 
by the triplet of values {Z(0) = i,B(0) = u,r(0) = ro}. To this end we can define the 
force of interest at the generic time s as follows: 

8i,u,r (s) = x(Ti > s\J = i,T = -u,Ti > 0)r i>ro (s) ^ 
+ x(Ti < s\J = i,T = -u,Ti > 0)5j lArijro ( Tl )(s - T x ). 

Here x(A | B) is the indicator function of set A given the information B. 

Remark 2.1. The processes x{T\ > s\Jo = i,Tq = —u,T\ > 0) and rj jro (s) are inde- 
pendent for any s € M+ and i G /. 

Remark 2.2. The process x(s < T\ < s + h\J$ = i,T^ = —u,T\ > 0)5j 1)0i r. 4 ro (Ti)( s + 
h — T\) conditioning to the value of Vj iro (Ti) is independent o/rj /ro (s). 

We are interested in the ZCB process Vi^ U)rt (t, T) := exp(— d~i, u ,r t { s )ds) expressing 
the value of a zero coupon bond paying 1 Euro at time T given that at current time 
t the semi-Markov process is in state i and it is entered in this state u periods before 
and the force of interest at that time is r^. To this end in what follow we would like 
to consider only diffusion processes that allow an explicit representation of the Laplace 
transform of f ri !ro (s)ds for example |Vasicek(1997)| , |Hull and White(1990)| and Cox, 
Ingersoll and Ross (1985) (CIR) model. 
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Theorem 2.3. Let 

V%] rt (t,T)=E[(v i , u , rt (t,T)r] (2.4) 
be the nth order moment of the ZCB process then it results that 

<t,(^) = 1 -f^ ) + " ) ( g W( t ,r)) 

where B^ is defined by 

B^ t (t,T) = E | exp (- n J% hrt ( S )ds^) } (2.6) 

Proof. Let consider, at the time t, the process Vi )U>rt (t,T) := exp(— Si )U ,r t { s )ds) and 
condition to the time of next transition of the switching process. We can partition the 
state space in two possible events {T N n\ +1 > T} or {T N n\ +1 < T}. 
The first event {Tjy(t)+i > T} corresponds to the possibility to have no transition up to 

the time T and it has probability 1 jz^j. ■ Under this event, the force of interest 
is given by 

Si, u ,n(s) = n,r t {s) for s G [t,T]. (2.7) 

In this case, the ZCB process assumes value 



v i,u,r t {pi T) 

T N(t)+l> T 

tacking now the expectation we get 



exp ^- r itn (s)dsj , (2.8) 



n(vi,u,r t (t,T)) n \T m+1 >T} = B™(t,T). (2.9) 
Then we have that 

E{(« i)Uirt (*,r)) n x(T JV(t)+1 >T\J N{t) =i,T N{t) =t-u,T m+1 >t,r(t) = r t )} 
- l - H ^- t + U) B^it,T). (2.10) 



1 - Hi(u) 



The other event {Tjv(i)+i G (t, T], Jjv(t)+i G ^} corresponds to the possibility to have 
at least one transition in the considered time interval. The probability that the semi- 
Markov process has the first transition in the time interval (r, r + dr), for r G (t,T), 
into a state k G /, is given by 

P(Mt)+i = k, T m+1 G (r, r + dr) | jf M , Z(t) = i, B(t) = u) = f^^ dr, 

where (J~t)t is the filtration generated by the process (Z, B), that is the semi-Markov 
and the backward recurrence time together. To properly evaluate the expectation under 
this event, we have to consider even all the possible value that the force of interest can 
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assume on transition time r. Denoting by F r . r ( T -j the cumulative distribution function 
associated to ri >rt (r), using the continuity property of the force of interest, we have that 

P{JN{t)+i = k,T N (t)+i G (r,r + dT),6 ijUtn (r) € (x , x + dx)\T t ) = JZ~jf(^\ F n,r t (r) (tfa)^T 

where 7^ is the filtration generated by (Z,B,W). Under the event ^4 = {Jjv(t)+i = 
k,Tff u\ + i € (r, r + dr), <5i )U)rt (r) G + dx), J^} the ZCB process can be expressed 
as 



v i,u,rt (ti T) 



exp 



0~i,u,rt (s)ds 



r iin (s)ds \exp - / <5fc,o,x(s)a!s 



Moreover from remark f2, 21 we have 



E{(w;, u ,rt(t,r)) n |yl} = E jexp I -n y r i)rt (s)da J exp i-n J 5 kfiiX (s)ds j \A 

= E^exp f-n J\ iirt (s)ds) {A^E^exp (-n £ 6 kt0>x {s)ds\ \a\ = B^(t,T)V^(r,T). 
and consequently 

E{(vi !U:n (t,T)) n x(t < T N{t)+1 <T\J N( t) =i,T N ( t ) = t-u,T N{t)+1 > t,r(t) = r t )} 



kei 



Hi{u) 



(2.11) 



The value of the zero coupon bond is given by 



E{(v i>u>rt (t,T)) n } = E{(v i>u>rt (t,T))^ 

+ E{(v i:U!rt (t,T)) n x(t < 2jv(t) + i < T\J N{t) = i,T N{t) = t - u,T N(t)+1 > t,r(t) = r t )}, 

(2.12) 

then by substitution of expressions (|2.10p and (I2.1ip in f|2. 12|) the proof is complete. 



□ 



Corollary 2.4. Let denote by Ri )U ,r t { s ) = ^[^i,u,r t { s )]> then it results that 



Ri,u,rt Ov 



l-Hj(s + u) 
l-Hi{u) 



here, m ijrt (s) := E[r itn (s)]. 



Rk,0,x( S ~ T ) F n,r t (r)( dx ) ) dT 

(2.13) 



We can derive the following result that help us on obtaining the covariance function 
of the force of interest. 



Theorem 2.5. Let s > and h > 0, and denote by 

Zi, u ,r {s,h) =E[6(s)5(s + h)\Z(0) =i,B(Q) =u,5(0) =r ]. 



(2.14) 
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Then it results that 

1 - HAs + ft + u 



i,u,r (s,h) — . Pi,ro{s,S + K) 



1-Hiiu) 

'■s+h 



, \- f S+n Qik(T + u) ^r°° R l ± h \W (A\A 

+ 2^ / l - H ( ) m *' r °( s ) / Rk,o,x{s + h - T)F iMT) (dx)dT ^ ^ 
+ / 1 ETTT" / s fc,oM s - t, s + ft - T)F iMr) (dx)dT 



kei 

here, p it r (s,s + h) := E[r itro (s)r^ TO (s + ft)]. 

Proof. Using the definition of force of interest process 8, cf. formula (|2.3p . we have 
that 

5i,u,r (s)5i :Utro (s + ft) = {x(Ti > s[J = *,T = -u,Ti > 0)r i)T . o (s) 

+ x(Ti < s |Jo = i 5 r = -u,r 1 >o)[5 JlAr . ro(Tl) ( s -r 1 )]} 

x { x (Ti > s + ft|J = i,T = -u,Ti > 0)rv (s + ft) 

+ < s + ft|J =i,T = -u,Ti > 0)[6j 1)0>ritro (T l )(s + h-T 1 )]}, 

solving the product we obtain three terms 

8i,u,ro(s)8i,u,r (s + h) = x(?i > s + ft|Jo =i,T = -u,Tx > 0)r ijro (s)r ijro (s + ft) 
+ X(« < Ti < s + h\J =i,T = -u,Tx > 0)r i)ro (s)5j li0iriro ( Tl )(s + h-T x ) 

+ x(Ti <s\J = i,T = -u,T t > 0)<5 JlAr<ro ( Tl )(s-T 1 )<5 JlAni , o ( Tl )(s + ft-T 1 ), 

(2-17) 

The first corresponds to the possibility that semi-Markov process has no transition up 
to time s + h, the second to the case to have the first transition into the time interval 
(s, s + ft) and the third takes into account for the possibility to have the first transition 
before s. 

Taking expectations of the first term on right hand side (r.h.s.) of (I2.17|) and nothing 
that (cf. remark [2.ip the random variable x(^i > s + ft| Jo = i, Tq = —u,T% > 0) is 
independent both from ri jrQ (s) and ri )TQ {s + ft) yields, 

E[x(Ti > s + ft|J = i,T = -u,Ti > 0)r i)ro (s)r ijrQ (s + ft)] 
= E[ X (Tx > s + ft| J = i,T = -u,Ti > 0)]E[r i>n) (a)r i)ro (« + ft)] (2 lg) 

l-Hi(s + h + u) . , 
1 - iJt( u ) 

Regarding the second term on r.h.s. of (|2.17p . tacking the expectation and con- 
ditioning on the value of rj /ro (Ti), we can use the independence between rj )7 . (s) and 
X(s < T i < s + ^|Jo =i,T = -u,Ti > 0)<5j lAriro ( Tl )(s + ft - Ti) (cf. remark to 
obtain 

E{E[ X (s < Ti < s + ft|J = i,T = -u,T x > 0)r hro (s)5 Jlfi ^ ro{Tl) (s + ft - Tx)\n, ro (Tx) = x]} 
= E{E[ X (T! > 8 + ft|J = i,T = -u,T x > 0)<5 JlAr ^ o(Tl) (s + h- T x )\r^ rQ {Tx) = x]}E[r i>ro (s)) 

= JT* ^^nn,r (s)} f£ E[S ki0 , x (s + ft - r)}F lMr) (dx)dT, 
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but 

E[r iiro (s)} = m ijro (s), (2.19) 

and 

E[S ki0 , x (s + h-r)}= R k , , x { s + h-r), (2.20) 
then, by substitution, we get 

E[ x (s <T l <s + h\J = i,T = -u,T x > 0)r ijro (s)S Jl>Oir .^ Tl) (s + h- Ti)] 

= Ekei I S s +h W$ m i,r (s) f£ R k , 0lX (s + h ~ r)F lMT) {dx)dr. (2.21) 
Finally, tacking the expectation of the third term on l.h.s. of (12, 17ft we obtain 
E[x(Ti < s\J = i,T = -u,Ti > 0)S Jl!O>r . ro (T 1 )(s - T^d^^.^^is + h- Ti)] 
= Efee/ Jo T=7^ E [ S kAri,r (T)(s ~ T)8 k>0 , nro(T) (s + h- r)]dr, 
and by considering all possible values of r- l)ro (r) we get 

E[x(Ti < s\J = i,T = -u,Ti > 0)6 Jl)0>r . ro ( Tl) (s - Ti)S Jl>0tri inj (n)(* + h ~ T i)} 
= Zkei So ( f£ m,°As - r) W + h - r)]F lMr) (dxj) dr, 

but 

H$k,o,x(s - T ) S k,o,x( s + h-r)] = E k! o >x (s -T,s + h-r), 

then 

E[x(Ti < s\J = i,T = -u,Ti > 0)<5 JlAriro(ri) (s - Ti)5 Jl0 |P . >rQ ( Tl )(s + h- 2\)] 
= E k ei So ¥01 ( S-™ S*^(- ~r,s + h- r)F tMT) (dx))dr, 

□ 

sectionParticular Models So far, we do not specified any specific dynamic driving the 
force of interest between two consecutive renewal moments. In this section we would like 
to consider some particular diffusion models for the force of interest as |Vasicek(19 97) |, 
[Hull and White(1990)| and |Cox et al.(1985)] . For these models, as well known, it is 
possible to obtain an explicit representation of the Laplace transform of the integral of 
the force of interest r. We will see as our general results specializes in these particular 
cases. Without loss of generality, as long as we are working with an homogeneous time 
model, we can suppose that the present time is t = 0, with Jo = i, then we will study 
the dynamic of the force of interest from up to a generic instant t < T\. 

2.1 Vasicek modulated model 

We assume that the process rj jro (i), satisfies 

dri, ro (t) = ai(bi - r ijTo (t))dt + aidW t 
»Vo(0) = r ° 
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where ctj, hi and <jj, for fixed i £ J, are non negative constants. The solution of the 
previous stochastic differential equation is 



»Vo(*) = &i + (n> - h)e °** + 



e a * h dW, 



(2.22) 



For any fixed t, ri iTQ (t) is normal distributed, ri >ro (t) ~ J\f(mi jro (t),af ro (t)) with mean 



m i} r (t) = E[r i)ro (i)] = 6j + (r - bi)e 



(2.23) 



and variance 



^ (*) = g:(i-e- 2c "' 



(2.24) 



It follows that (see Lamberton and Lapeyre(2008)| ) Jq r(s)ds is a normal random vari- 
able. Then we can completely characterize the distribution of J Q ri^ ro (s)ds with its mean 
and variance. The mean, using Fubini's Theorem, is given by 



E 



o 



r it r (s)ds 



E[r itro (s)]ds = I m i , ro (a)cfa = 6it + ^_Ji(i_ e -«<t) j (2.25) 



and the variance by 



Var 



/ r iiro (s)ds 
Jo 



a i a i 



(1-e 



-ait\ 



-W(l-e 



(2.26) 



Now, we can give an analytical expression of the n moment of a ZCB value between 
two consecutive renewal moment, B^} (0, t) (cf. theorem I2.3p . defined by 



B<J>(0,t)=E 



exp{— n / rj jf , (t)dt} 



(2.27) 



which can be seen as the Laplace transform of the Normal random variable ri^ ro (t)dt 
and then it can be expressed as 



n t r (t)dt 



+ Y Var 



r i<ro (t)dt 



(2.28) 



Substituting the expression for mean and variance we get 

'ofra 2 n(r - h) 



B^ Q (0,t) = ex P 



afn 2 



nbi I t 



+ 



cr?n 2 



(l-e-^)-^(l-e 



-aA \2 



2af 



Now, we are able to express all the moments of the ZCB in any time interval. 



(2.29) 



Corollary 2.6. The nth order moment of the ZCB process described by a Vasicek model 
modulated by a semi-Markov process of kernel Q is 



»5<S,(«,T)-EW-n / .!,„„(«)•'■«) 



T 



l-H t (T-t + u) (n) 

l-Hiiu) ^WMJJ 



+ E [ (*, r) £j f mm , n ^ t)<n (r _ t)) (x)^ (r, T)d.) dr. 

(2.30) 
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where B^(t,T) is given in equation \2. 30\) and /^( mj r ( T -t),<r? (r-t)) * s ^ e distribution 
function of a normal with mean nii >rt (r — t) an<i variance af (t — t) ( cf. eq. \2.2'J\) and 
\2. 24\ ) respectively). 

Corollary 2.7. The first moment of the force of interest of the Vasicek model modulated 
by a semi-Markov process of kernel Q is 



(2.31) 



where fj^f( m . r ( T ) CT 2 (t)) * s ^ e distribution function of a normal with mean mi >rt (r) and 
variance erf (t) (cf. eq. h2.23\) and \2. 24\ ) respectively) . 

Corollary 2.8. The product moment of the force of interest of a Vasicek model modu- 
lated by a semi-Markov process of kernel Q is 

Z iiU , ro (s, h) = E[6(s)6(s + h)\Z(0) = i, B(0) = u, 6(0) = r ] = ^—^±^±A pi ( s , s + h) 
+ E / T- h J + ^ - / ^K„W^, ro W)W^A^ + h - r)dxdr 

+ E / TT^Y / &/"(«K,n> (r),o?, ni (r))^)Sfc I o,x( s t } s -\- h r)dxdr 

(2.32) 

where Pi t r (s, s + h) is given by 

Pi, ro {s,s + h) = E{n, ro (s)r hro (s + fc)] = <rfe~ a< ' 1 , (2.33) 

and //\^( mirt ( r ) o-2 ( T )) is the distribution function of a normal with mean rn,j in (r) and 
variance o~f (t) (cf. eq. \2.23\) and {2. 24\ ) respectively) . 

We conclude this subsection with a Monte Carlo algorithm able to generate the 
trajectories of the Vasicek model modulated by a semi-Markov process of kernel Q in 
the time interval [0, T]. The algorithm consists in repeated random sampling to compute 
successive visited states of the random variables {Jo, Ji, ...}, the jump times {Tq,Ti, ...} 
and the force of interest process r(t) up to the time T. 
The algorithm consists of 5 steps: 

1. Set n = 0, Jo = i, Tq = 0, r(0) = ro, horizon time= T; discretization step= h 

2. Sample J from pj nt . and set J n +i = J(u); 

3. Sample W from Gj nt j n+l (-) and set T n+ i =T n + W(uj); 

4. For each i = T n : h : max(T, X)u +1 — 1) 
Sample TV from A/"(0, h) 

Set r(i + 1) = b Jn + e~ a ^ (r(i) - bjj + e~ aj "a Jn N(i + 1) 

5. if T n +i > T stop 

else Set n = n + 1 and go to 2). 
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2.2 Hull and White modulated model 

Let's now assume that the process rj >ro (i) satisfies 

dn, ro {t) = (a f (t) - Pi(t)r ijro (t))dt + Gi{t)dW t 
^,r (0) = r 

where aj, /3j and cij are deterministic function of time. The solution of this stochastic 
differential equation can be expressed (see |Shreve(2004)| ) as 



Ti,ro (^) — ^ 



-ki(t) 



r + e k ^a i ( y u)du + / e k ^a t (u)dW(u) 



where 



ki(t) = / f3i(u)du. 
Jo 

Note that (?"i,ro(^))te[0,T] i s a Gaussian process whose mean is given by 

mi,r (t) = e~ h{t) r + [ e k ^ u) ai{u)du , 

J o 

and its variance is given by 

( t ) =e - 2 *iW f e^a^ufdu. 
Jo 

Morover the process (Jq ?"i,ro( s )^ s )te[o,T] is Gaussian as well, with mean 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



E 



r it r (s)ds 



-ki(t) 



ro+ e fcl(u) ai(u)du 



ds, 



and variance 



Var / r iro (a)ds 



2ki{u)J2 



*/(«) 



- k ^ds\ du. 



(2.38) 



(2.39) 



Then, we are able to express the n th moment of a ZCB value between two consecutive 
renewal moment in an analytical form 



Bjj>(0,t)=E 



exp{—n / ri nt {t)dt} 



exp < — E 



r 



r i<ro (t)dt 



+ -Var 



n tro (t)dt 



(2.40) 



In analogy to the Vasicek model we could give similar results for the value of the n th 
moment of ZCB in an arbitrary time interval and the first two moments of the force of 
interest. We omit here the details for the sake of brevity. 
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2.3 CIR modulated model 

In this model we assume that the dynamics of the force of interest is described, be- 
tween two consecutive renewal moments, by the following stochastic differential equation 



dri, ro {t) = (oj - bir it r (t))dt + ai^/r~Jt)dW t 
r i,r (0) = r o 

here, cij and <jj are non negative constants while hi 6 R. This stochastic differential 
equation, for any ro £ R, has a unique solution (see Ikeda and Watanabe(1988)| ). 



We will not write the explicit solution of this equation, but following the approach 
of |Lamberton and Lapeyre(2008)| , we will study the property of the distribution of 



r 



and 



;,r (i)- To this aim, let first study the distribution of {ri iTo (t), f r{ t r Q (s)ds). It is 
possible to show (see |Lamberton and Lapeyre(2008)] ) that 

E (e-ArW^e-M/o'W^) = exp(-a^, A ,^))exp(-r VwW) (2-41) 
where the functions <pi.\^{t) and i/Ji t \ t u{t) are given by 

2 ( 2 e (l ^ hl) \ 

erf ^of A(eT** - 1) + 7i - 6j + e^*(7j + 6j) y 

, m _ A( 7 i + bi + e^(7i-bi)) + 2//(e^-l) 
M " o?A(tf»* - 1) + 7 i - 6| + + k) [ZA6) 

with 7j = yjbf + 2af fi . Using this result we can explicitly evaluate the Laplace trans- 
form of Ti ro (t), indeed by putting fj, = in the above expression we get 

Moreover, we can obtain an analytic representation of the n*' 1 moment of a ZCB value 
between two renewal moments, in fact setting A = and fi = n we have 

( 2 7i e i ^ \^ / 2(e**-l) \ 
- exp -r \ • rr ) ] (2.45) 



7i - 6i + e^*(7i + 6i) y V H-h + e^^i + bi 



where 7j = ybf + 2<rfn . Starting from here we can obtain similar result as for the 

Vasicek model for the n th moment of a ZCB in an arbitrary time interval and the first 
two moments of the force of interest. We omit here the details for the sake of brevity. 



13 



References 

[Bj6rk(2004)] Bjork, T., Arbitrage Theory in Continuous Time, Oxford University Press 
Inc, New York (2004). 

[Cox et al.(1985)] Cox, J.C., Ingersoll, J.E., and Ross, S.A., A theory of term structure 
of interest rates. Econometrica. 53, pp. 385-407 (1985). 

[D'Amico et al.(2009)] G. D'Amico, J. Janssen, and R. Manca, Semi-Markov Reliability 
Models with Recurrence Times and Credit Rating Applications, Journal of Applied 
Mathematics and Decision Sciences, Article ID 625712, 17 pages (2009). 

[D'Amico et al.(2011)] G. D'Amico, J. Janssen, and R. Manca, Duration Dependent 
Semi-Markov Models, Applied Mathematical Sciences, Vol. 5, no. 42, 2097 - 2108 
(2011). 

[Duffle and Kan(1996)] Duffie, D., Kan, R., A yield factor model of interest rates. Math. 
Finance, 5, 379-406 (1996). 

[Hull and White(1990)] Hull, J., White, A., Pricing Interest Rate Derivative Securities. 
The Review of Financial Studies, vol. 3, n. 4, pp. 573-592 (1990). 

[Hunt and Devolder(2011)] Hunt, J., Devolder, P., Semi-Markov regime switching in- 
terest rate models and minimal entropy measure. Physica A, 390(21-22), 3767-3781 
(2011). 

[Ikeda and Watanabe(1988)] N. Ikcda and S. Watanabe. Stochastic Differential Equa- 
tions and Diffusion Process. North-Holland (1988). 

[Janssen and Manca(2007)] J. Janssen and R. Manca, Semi-Markov risk models for 
Finance, Insurance and Reliability, Springer, New York, (2007). 

[Lamberton and Lapeyre(2008)] D. Lamberton and B. Lapeyre. Introduction to 
Stochastic Calculus Applied to Finance. Chapman & Hall/CRC (2008). 

[Mamon(2002)] Mamon, R.S., "A time-varying Markov Chain model of term structure" 
Statistics and Probability Letters, 60 pp. 309-312 (2002). 

[Mamon(2004)] Mamon, R.S., "Analytic pricing solutions to term structure derivatives 
in a Markov chain market" IMA Journal of Management Mathematics, 15, pp. 
243-252 (2004). 

[Norberg(1995)] Norberg, R., "A Time-Continuos Markov Chain Interest Model with 
Applications to Insurance." Applied Stochastic Models and Data Analysis, 11, pp. 
245-256 (1995). 

[Pye(1966)] Pye, C, "A Markov model of term structure." Quarterly Journal of Eco- 
nomics, 25, pp. 60-72 (1966). 

[Rebonato et al.(2005)] Rebonato, R., Mahal, S., Joshi, M., Buchholz, L.D., Nyholm, 
K., Evolving Yield Curves in the real- World Measures: A Semi-Parametric Ap- 
proach. J. Risk 7(3), 29-62 (2005). 



14 



[Shreve(2004)] S. Shreve. Stochastic Calculus for Finance II: Continuous-Time Models. 
Springer (2004). 

[Stenberg et al.(2007)] Stenberg, F., Manca, R., Silvestrov D., An algorithmic approach 
to discrete time non-homogeneous backward semi-Markov reward processes with an 
application to disability insurance. Methodol. Comput. Appl. Probab. 9, 497-519 
(2007). 



[Vasicek(1997)] Vasicek, O., "An Equilibrium Characterization of the Term structure." 
Journal of Financial Economics, 5, pp. 177-188 (1977). 



